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ABSTRACT 



We have discovered both intermediate-order gravity mode and low-order pressure 
mode pulsation in the same star, HD 209295. It is therefore both a 7 Doradus and a 
S Scuti star, which makes it the first pulsating star to be a member of two classes. 

The analysis of our 128 h of multi-site spectroscopic observations carried out 
over two seasons reveals that the star is a single-lined spectroscopic binary with an 
orbital period of 3.10575 ± 0.00010 d and an eccentricity of 0.352 ± 0.011. Only weak 
pulsational signals arc found in both the radial velocity and line-profile variations, but 
we have succeeded in showing that the two highest-amplitude 7 Doradus pulsation 
modes are consistent with 1=1 and \m\ = 1. 

These two modes dominated our 280 h of BVI C multi-site photometry, also ob- 
tained over two seasons. We detected altogether ten frequencies in the light variations, 
one in the 5 Scuti regime and nine in the 7 Doradus domain. Five of the 7 Doradus 
frequencies are exact integer multiples of the orbital frequency. This observation leads 
us to suspect they are tidally excited. Attempts to identify modes from the multicolour 
photometry failed. 

We performed model calculations and a stability analysis of the pulsations. The 
frequency range in which 6 Scuti modes are excited agree well with observations. 
However, our models do not show excitation of 7 Doradus pulsations, although the 
damping is smaller in the observed range. We also investigated tidal excitation of 
7 Doradus modes. Some of the observed harmonics of the orbital period were found 
to be unstable. The observed orbital harmonics which arc stable in the models can be 
understood as linear combinations of the unstable modes. 

We could not detect the secondary component of the system in infrared photom- 
etry, suggesting that it may not be a main-sequence star. Archival data of this star 
shows that it has a strong ultraviolet excess, the origin of which is not known. The 
orbit of the primary is consistent with a secondary mass of M > 1.04Mq, which is 
indicative of a neutron star, although a white dwarf companion is not ruled out. 

Key words: stars: variables: 6 Set - stars: oscillations - stars: individual: HD 209295 
- stars: binaries: close - stars: binaries: spectroscopic - stars: neutron 
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1 INTRODUCTION 

Four different classes of multi-mode pulsating variables are 
found near the intersection of the classical instability strip 
and the main sequence. However, asteroseismology (probing 
the stellar interior through such pulsations) of these stars 
has proven to be more difficult than expected. Nonetheless, 
the example of helioseismology and its enormous reward in 
terms of the physical understanding of the Sun's interior 
(e.g. see Gough 2000) has been a great motivation for con- 
tinuing efforts in probing the structure of these stars by 
means of their pulsational properties. 

Many attempts to detect solar-type oscillations in stars 
have resulted in inconclusive results, although success has 
been claimed in some instances (e.g. Kjeldsen et al. 1995). 
Because of their extremely low amplitude, it is only recently 
(e.g. Bedding et al. 2002) that we are beginning to see detec- 
tions which may be considered significant. Successful aster- 
oseismology based on the analysis of solar-type oscillations 
now appears immanent. 

The rapidly oscillating Ap stars (see Kurtz & Martinez 
2000 for a recent comprehensive review) are high-radial or- 
der p-mode pulsators. Their pulsation spectra resemble that 
of the Sun. In these stars, only a few independent modes are 
observed. Moreover, the presence of a strong magnetic field 
is expected to modify the frequencies and eigenfunctions of 
the modes. Successful modelling of these pulsations does not 
seem possible at this time. Nevertheless, some progress is 
currently being made (e.g. Cunha & Gough 2000) . 

The third group of pulsators in this region of the HR Di- 
agram are the 5 Scuti stars, which are low-radial order p- and 
probably mixed-mode pulsators (see Breger & Montgomery 
2000). For some of these stars, tens of pulsations modes 
have been detected and their frequencies determined with 
high precision (see e.g. Handler et al. 2000). Although mean 
densities, and even asteroseismological distances (Handler 
et al. 1997) have been estimated, the problem of mode iden- 
tification and shortcomings in the stellar models are major 
difficulties (see e.g. Pamyatnykh et al. 1998 for a discussion) . 

The 7 Doradus stars are high-overtone gravity (g) mode 
pulsators (Kaye et al. 1999a). In these stars only a few modes 
are excited to observable amplitude amongst a dense forest 
of possible modes. As a result, accurate starting values of 
the basic stellar parameters are required to enable observed 
and calculated frequencies to be matched. 

In principle, it is possible for 5 Scuti and 7 Doradus 
pulsations to co-exist in a star as they occupy overlapping 
regions in the HR diagram. This offers new opportunities for 
successful asteroseismology, as noted by Handler (1999a). In 
such stars it may be possible to use the 8 Scuti pulsations 
to place strong constraints on the stellar parameters, easing 
the mode identification problem of the 7 Dor pulsations. In 
this way it may be possible to obtain information about the 
deep interior which is unobtainable for pure 5 Set stars. 

Consequently, Handler & Shobbrook (2002) searched for 
S Scuti pulsations in all known candidate 7 Doradus stars 
located within the 5 Scuti instability strip and accessible 
from intermediate southern geographical latitudes. About 
one third of all non-Am and non-Ap stars in the lower in- 
stability strip are indeed 5 Scuti stars (Breger 1975). Be- 
cause Am and Ap stars are rare (if not absent) amongst the 
7 Doradus stars (Handler 1999a), detection of some "hybrid" 



Table 1. Multi-frequency solution for the 7 Doradus variabil- 
ity of HD 209295 from Hipparcos observations. Formal error esti- 
mates were derived from Montgomery & O'Donoghuc (1999). The 
signal-to-noise ratio was calculated following Breger et al. (1993); 
a S/N > 4 corresponds to a significant detection. 



ID 


Frequency 


Hp Amplitude 


S/N 




(cycles / day) 


(mmag) 




fi 


1.12957 ± 0.00002 


37 ± 2 


12 


h 


2.30222 ± 0.00003 


31 ± 2 


10 


h 


2.57579 ± 0.00005 


16 ± 2 


5 


f* = h - fi 


1.1726 ± 0.0001 


12 ± 2 


4 



Table 2. Published uvby/3 photometry of HD 209295. The differ- 
ences in the magnitudes and photometric indices are consistent 
with stellar pulsation. When brighter the star it is bluer and the 
ci and indices are larger. 



V 


b-y 


mi 


Cl 


P 


Reference 


7.32 
7.29 


0.149 
0.139 


0.182 
0.185 


0.819 
0.840 


2.781 
2.821 


Twarog (1980) 
Handler (1999b) 



stars may be expected if 5 Scuti and 7 Doradus pulsations 
are not mutually exclusive. In this paper we report the dis- 
covery and analysis of such a star, HD 209295 = Sangam 
Mani. 

This is a V=7.3 mag Southern (<5 = —64°) star of spec- 
tral type A9/F0 V (Houk & Cowley 1975). It was discovered 
as variable by the Hipparcos mission (ESA 1997). Handler 
(1999a) performed a frequency analysis of these data and de- 
tected four periods; we show the results in Fig. 1 and Table 
1. We are aware that a search for multiple periods in Hippar- 
cos photometry is difficult and can lead to spurious results 
(Eyer & Grenon 1998). However, the case of HD 209295 is so 
simple and convincing that an error in the analysis is con- 
sidered improbable. The multi-periodicity leaves no doubt 
that HD 209295 is a bona fide 7 Doradus star. 

HD 209295 is considerably hotter than all other 
7 Dor stars. Handler (1999a) suggested that the published 
Stromgren {b — y) of this star (Twarog 1980) is too blue com- 
pared to (B — V). Handler (1999b) obtained new Stromgren 
photometry and confirmed the discrepancy (cf. Table 2). 
Grenier et al. (1999) made two measurements of the radial 
velocity of the star using an echelle spectrograph and the 
cross-correlation technique. They did not comment on bina- 
rity or any other peculiarity. 

Using the mean of the Stromgren indices from Table 
2 (which indicate zero reddening), the effective temperature 
and gravity of HD 209295 can be estimated. Kurucz's (1991) 
calibration yields T eS = 7750 ± 100 K, log 3 = 4.10 ± 0.05. 

The Hipparcos parallax of HD 209295 is 8. 19 ±0.72 mas. 
Using Vq = 7.32±0.01, BC = -0.01±0.01, the parallax gives 
log L/L = 1.15±0.08. With logT off = 3.889 as determined 
above, R/Rq = 2.08 ± 0.09 and, by means of evolutionary 
tracks calculated using the Warsaw-New Jersey code (see 
e.g. Pamyatnykh et al. 1998), a mass of 1.84 ±O.O7M is 
determined. We can then also calculate a semi-independent 
second value for \ogg, namely 4.07 ± 0.05, consistent with 
the previous value. HD 209295 is therefore the hottest 7 Do- 
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Figure 1. Spectral window and amplitude spectra of Hipparcos 
photometry of HD 209295 with consecutive prcwhitcning of de- 
tected frequencies. The combination peak f2-fi is only considered 
significant because it occurs at a "predicted" frequency. 

radus star known to date; it lies in the middle of the 8 Scuti 
instability strip on the main sequence (Fig. 2). 

Because of its location in the colour-magnitude dia- 
gram, HD 209295 is a good candidate 6 Scuti star. Indeed, 
multi-periodic 5 Scuti pulsations were found by Handler & 
Shobbrook (2002), thus placing HD 209295 as a member of 
two classes of pulsating star. Spectroscopic and photomet- 
ric follow-up observations were organized immediately after 
this discovery, and a multi-site campaign was conducted in 
August 2000. We report the results of these campaigns. 



2 OBSERVATIONS 
2.1 Optical Photometry 

As comparison stars, we used HD 207802 (B9/B9.5V, V = 
8.0) and HD 209829 (F3V, V = 7.7) in our photoelec- 
tric photometry. Measurements were obtained at the 1.0-m, 
0.75-m and 0.5-m telescopes of the Sutherland station of 
the South African Astronomical Observatory (SAAO), the 
0.6-m telescope at Siding Spring Observatory (SSO) in Aus- 
tralia (Shobbrook 2000), and the 0.6-m Zeiss telescope at 
the Osservatorio do Pico dos Dias (LNA, Brazil). Observa- 




0.1 0.2 0.3 

(b-y) or (b-y) 
Figure 2. The location of HD 209295 (star symbol) in the colour- 
magnitude diagram compared to 252 known <5 Set stars (open 
circles) with available Stromgren photometry from the catalogue 
of Rodriguez, Lopez-Gonzalez & Lopez de Coca (2000). The blue 
and red edges of the 8 Set instability strip (Rodriguez & Breger 
2001) arc indicated. 



tions were conducted from October 1999 to October 2000. 
Most measurements were obtained in August and Septem- 
ber 2000 in a co-ordinated multi-site effort. A journal of the 
observations is given in Table 3. 

We used the Johnson B and V filters as well as the 
Cousins I c filter (the latter only in the year 2000) with a 
total integration time of ~1 minute in each filter as a com- 
promise between good time resolution and maximum colour 
information. Apertures of 30 - 45 arcseconds on the sky were 
used. Sky measurements were taken at suitable intervals de- 
pending on the brightness and proximity of the Moon. 

The observing sequence was chosen to result in both 
good coverage for the short-period S Scuti pulsations and 
best long-term stability. We adopted the sequence C1-V-C2- 
V-Cl... (CI and C2 are the comparison stars and the V is the 
variable), which resulted in one variable star measurement 
every ~6 minutes. Supplementary observations of UBV(RI) C 
standard stars were also acquired. 

Data reduction was performed in the standard way, i.e. 
corrections for coincidence losses, sky background and ex- 
tinction were followed by calculating differential magnitudes 
between the comparison stars. The latter were examined for 
variability. It seems that one of the two stars is slightly vari- 
able on a time scale of about 1.1 or 10 days with a V am- 
plitude of 2 mmag. We investigated whether this could be 
due to differential colour extinction and found no support 
for such an interpretation. As our programme star exhibits 
long-term variability with much higher amplitude, we can- 
not determine unambiguously which of the two stars is the 
potential variable, but we suspect HD 209829. 

In any case, we proceeded by constructing a differential 
target star light curve relative to the measurements of both 
comparison stars. We standardized these magnitude differ- 
ences by using the transformation equation slopes calcu- 
lated from the standard star observations mentioned above. 
We note that the photometric zero-points of the different 
telescope/instrument combinations agreed to better than 2 
mmag for each filter or colour used during the multi-site 
campaign, but we experienced some difficulties with homog- 
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Table 3. Journal of the photometric observations. The measure- 
ments before HJD 2451600 were made in the Johnson BV system, 
whereas later observations utilized Johnson-Cousins BVI C except 
for one night of BV observations marked with an asterisk. 



Telescope 


Run start 


Run length 


Observer 




HJD - 2450000 


(h) 




S A AO fl 7^ m 
Oii.lU U. 1 O -III 


1 4fi4 


3.31 


GH 


cja AD 7^ m 

On/iv U. I O -111 


1 4fi^ 94Q 


6.34 


GH 


QQf) fi m 
OOW U.U III 


1 4fi7 0^3 
11U 1 . uoo 


2.78 


RRS 


SSO fi m 

OOW U.U III 


1 4fik Q1 fi 


6.05 


RRS 


SSO 0.6-m 


1491.937 


2.35 


RRS 


SSO 0.6-m 


1498.918 


2.35 


RRS 


SAAO 0.5-m 


1502.282 


1.85 


GH 


SAAO 0.5-m 


1503.268 


2.62 


GH 


SAAO 0.5-m 


1505.267 


2.57 


GH 


SAAO 0.5-m 


1507.310 


1.08 


GH 


SAAO 0.5-m 


1509.295 


1.63 


GH 


SSO 0.6-m 


1510.924 


1.85 


RRS 


SAAO 0.5-m 


1511.264 


2.40 


GH 


SSO 0.6-m 


1512.933 


2.18 


RRS 


SSO 0.6-m 


1513.926 


0.94 


RRS 


SSO 0.6-m 


1516.930 


2.04 


RRS 


SSO 0.6-m 


1517.921 


1.20 


RRS 


SAAO 0.5-m 


1518.266 


0.86 


LE 


SAAO 0.5-m 


1519.283 


0.48 


LE 


SSO 0.6-m 


1524.930 


1.63 


RRS 


«qn fi m 

DOU U.U III 


1754.110 


3.41 


RRS 


SSO fi-m 

k_> kj v_y u.u in 


1 754 962 


5.45 


RRS 


SAAO fl 7^ m 

Orv/lU U. 1 O -111 


1 7^.8 343 

.1 f UO.OIO 


3.46 


ERC 


SSO fi-TYl 

OUU U.U 111 


1759.947 


8.11 


RRS 


SAAO (1 7^ m 

On/iU U. I U 111 


17fi1 343 

1 1 Ul. OIO 


7.87 


ERC 


SAAO 7^-™ 

ijjinu u. i u in 


1 7fi2 3fi4 

JL I U__- . -JUrt 


7.25 


ERC 


SSO fi in 

DOW U.U III 


1 7fi9 Q44 


7.75 


RRS 


S A AO 75-m 

ijjirvu u. i u in 


1 763 337 

JL 1 U -J. <-)•) 1 


8.11 


ERC 


S A AO fl 7^ m 

O/irtU U. I U III 


1 7fi4 978 


6.53 


ERC 


SSO fi-m 

k_> kj V_/ U • U 111 


1 7fi7 227 

1 1 U 1 .i-ii 1 


2.28 


RRS 


S A AO n 7^ m 

Jrt/lvJ U. i J 111 


1 7fi8 941 

JL ( UO. i4 JL 


10.25 


GH 


SAAO 7^ rn 
O/i.lv U. 1 O -111 


1 7fiQ 9^Q 


10.27 


GH 


l_jlN.fl. U.U 111 


1 7fiQ 81 Q 


0.24 


AB 


SAAO 7^ m 

OnHU U. 1 O 111 


1 770 9^4 


10.37 


GH 


T,NA fi m 

IjIn _ri U.U 111 


1770.514 


6.84 


AB 


S A AO D 7^ m 

O/i.lU U. 1 U 111 


1771.237 


10.27 


GH 


T,NA fi-m 

lJlN.rt U.U 111 


1771.498 


3.19 


AB 


S A AO f) 7^ m 

On/iv U. 1 O 111 


1 779 9^8 

.1 f I Zi.-jOO 


10.27 


GH 


S A AO D 7^ m 

Oil-lU U. 1 U 111 


1 77^ 99fi 


10.49 


GH 


SAAO 7^ m 

O/i.lU U. 1 U 111 


1 77^ 993 


4.99 


GH 


TWA fi m 

JLjlN.fl. U.U 111 


1 77^ ^7^ 


5.59 


AB 


TWA fi m 
l_.IN.rl. U.U 111 


1 77fi fi^l 

1. / 1 U.UiJl 


3.91 


AB 


TWA fi m 

l__TN.fl. U.U 111 


1 777 48fi 


7.56 


AB 


SSO D fi m 

OOW U.U III 


1 777 Q04 


6.14 


RRS 


OAAU U. (O-Ill 


1 778 4Q^ 


3 fi7 


\jrl 


SAAO 1.0-m* 


1785.256 


9.94 


GH 


SSO 0.6-m 


1785.918 


8.45 


RRS 


SSO 0.6-m 


1786.991 


0.67 


RRS 


SSO 0.6-m 


1787.898 


9.58 


RRS 


SAAO 0.5-m 


1789.242 


8.74 


GH 


SSO 0.6-m 


1789.935 


8.52 


RRS 


SAAO 0.5-m 


1791.337 


1.27 


GH 


SAAO 0.5-m 


1792.253 


8.42 


GH 


SAAO 0.5-m 


1850.284 


1.51 


DJJ 


SAAO 0.5-m 


1852.277 


1.78 


DJJ 


SAAO 0.5-m 


1853.332 


0.41 


DJJ 


SAAO 0.5-m 


1857.270 


2.45 


DJJ 


SAAO 0.5-m 


1858.276 


1.82 


DJJ 


SAAO 0.5-m 


1859.274 


1.92 


DJJ 


SAAO 0.5-m 


1862.278 


2.04 


DJJ 



Table 4. Standard optical and infrared magnitudes of HD 209295 
as measured on HJD 2451792.462. 



Filter 


Magnitude 


B 


7.623 ± 0.005 


V 


7.370 ± 0.005 


Ic 


7.070 ± 0.005 


J 


6.872 ± 0.009 


H 


6.727 ± 0.007 


K 


6.746 ± 0.008 


L 


6.721 ± 0.081 



enizing the data from the discovery season. Finally, the time 
base of our observations was converted to Heliocentric Julian 
Date (HJD). Most of the light curves obtained during the 
central part of the multi-site campaign are shown in Fig. 3. 

2.2 Infrared Photometry 

In addition to the optical data, L. A. Crause obtained in- 
frared JHKL measurements on the night of September 1, 
2000 using the 0.75-m telescope of the SAAO with the Mkll 
infrared photometer (an upgraded version of the instrument 
described by Glass 1973). The data were reduced to the 
SAAO system by using standard stars defined by Carter 
(1990). As the infrared observations were obtained simulta- 
neously with optical photometry, standard BVI C JHKL mag- 
nitudes could be calculated (Table 4). The star was ap- 
proaching a local minimum in its light curve when these 
measurements were taken. 

These optical and infrared colours are consistent with 
the spectral type of the star as inferred from standard rela- 
tions (Drilling & Landolt 2000, Tokunaga 2000), but agree 
less well with the effective temperature determined from 
uvby/3 photometry. However, this may be due to the high 
amplitude of the light variations. The total colour amplitude 
implies temperature variations of ?a 500 K between light ex- 
trema, sufficiently large to explain the apparent discrepancy. 

2.3 Spectroscopic observations 

Spectroscopic observations of HD 209295 were obtained with 
three different telescopes. We used the 1.9-m telescope at 
SAAO, the 1.9-m telescope at Mt. Stromlo (MS) in Aus- 
tralia, and the 1.2-m Euler telescope at the European South- 
ern Observatory (ESO) in Chile. A summary of these mea- 
surements can be found in Table 5. 

The bulk of the data originated from SAAO using 
the GIRAFFE echelle fibre-fed spectrograph attached to the 
Cassegrain focus of the 1.9-m telescope. The GIRAFFE spec- 
trograph has a resolving power of about 32 000. The 1024 x 
1024 TEK CCD chip gives a resolution of 0.06 - 0.09 A per 
pixel. A Th-Ar arc lamp was used for wavelength calibration 
with arc spectra taken at regular intervals to calibrate pos- 
sible drifts. The wavelength range was 4400 - 6680 A spread 
over 45 orders. Exposure times were normally 10 min for 
a S/N ratio of about 30 - 60. A total of 422 spectra of 
HD 209295 was obtained during three observing runs. 

The spectroscopic observations at MS were obtained 
with the coude echelle spectrograph on the Mount Stromlo 



Simultaneous p-mode and g-mode pulsation of HD 209295 5 



-0.4 
-0.3 
-0.2 

-0.4 
-0.3 
-0.2 

-0.4 
-0.3 
-0.2 

-0.4 
-0.3 
-0.2 

-0.4 
-0.3 
-0.2 

-0.4 
-0.3 
-0.2 

-0.4 
-0.3 
-0.2 

-0.4 
-0.3 
-0.2 




70.3 70.4 70.5 70.6 "^70. 7 70.8 71.2 71.3 71.4"" 71.5 71.6 



72.3 72.4 "^72.5 72.6 73.2 73.3 



....... •• ... -M 

....... *. ... *••» 

■ ■■■■■■ 




••*...•...•. 



76.7 76.8 77.5 77.6 -^77. 7 77.8 77.9 



_L 



73.4 73.5 



i .... i . 



73.6 75.5" 



78 



75.6 75.7 




78.1 



78.5 78.6 



_L 



■ I .... I 



85.3 85.4 85.5 85.6 



86 86.1 86.2 
HJD 2451700 + 



9.3 89.4 89.5 89.6 



Figure 3. Light curves of HD 209295 obtained during the multi-site campaign in August 2000. The asterisks are measurements from 
SSO, the open circles those from SAAO and the filled circles are LNA data. The upper light curves in each panel are the V data, the 
(V — I c ) variations are shown in the middle and the lowest are the (B — V) light curves. Zero-points are relative to HD 209829; the V 
zero-point is shifted by +0.11 mag for a better display. Note the multi-periodic slow and rapid light variability and the corresponding 
colour changes. 
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1.9-m telescope. Spectra with a resolution of ~ 0.15 A at Ha 
were recorded on a 2048x4096 pixel TEK CCD. The wave- 
length range obtained with this set-up was 4500 - 6900 A 
spread over 42 orders. Exposure times were between 17.5-20 
min for a S/N ratio of about 25-60. A Th-Ar arc lamp was 
used for wavelength calibration. Due to poor weather a total 
of only 15 spectra of HD 209295 was obtained. 

The spectroscopic observations in Chile were performed 
with the Swiss 1.2-m Ritchey-Chretien Euler telescope at 
ESO, La Silla. Euler is equipped with a high-resolution 
echelle spectrograph, CORALIE, and a 2k x 2k CCD camera 
with 15pm pixels. The resolving power amounts to 50 000 
and the total wavelength range is 3900A - 6800A in 68 or- 
ders, without any gaps in the coverage. The CORALIE spectra 
are extracted on-line following a standard echelle reduction 
procedure. In the case of our measurements, wavelength cal- 
ibration utilized the most recently obtained ThAr spectrum. 
Each order of the stellar spectrum was then divided by the 
blaze function. For a full description of the reduction scheme 
we refer to Baranne et al. (1996). The integration times were 
adapted to the atmospheric conditions (seeing, presence of 
clouds) and ranged from about 15 to 18 minutes. The typi- 
cal S/N ratio was « 30. In total, 38 spectra were obtained 
during one week. 

In contrast to the photometric observations, the spec- 
troscopic measurements were not reduced centrally. For 
the SAAO data, we used the local red uction software (see 
http://www.saao.ac.za/facilities/) for corrections for 



Table 5. Log of the spectroscopic observations of HD 209295. 
The date with respect to HJD 2450000, the run lengths, the mean 
integration times t, and the number of spectra, TV, are given. 



for bias and flat field, order extraction and wavelength cal- 
ibration. The same reduction steps for the MS data were 
performed in IRAF. 



3 ANALYSIS 

3.1 The photometry 

Our frequency analysis was performed with the pro- 
gramme PERIOD 98 (Sperl 1998). This package applies 
single-frequency power spectrum analysis and simultaneous 
mult i- frequency sine-wave fitting, but also has some ad- 
vanced options which will be described later as necessary. 
We calculated the spectral window and amplitude spectra 
of our data as well as the amplitude spectra of residual light 
curves after the previously identified periodicities had been 
removed using a multi-periodic fitting algorithm. Similar 
analysis were performed for all the three filters used. We 
adopted the mean frequencies for our final solution. 

Our analysis was concentrated on multi-site data from 
the year 2000. Due to the different telescopes and detectors 
used, the 1999 data analysis is complicated by zero-point 
calibration problems. Also, the data runs were rather short 
in 1999 and the possibility exists that the star may have had 
a different amplitude at this time. We used the 1999 data 
only when the zero-point problems were relatively minor. 



3.1.1 The 7 Doradus pulsations 

In Fig. 4 we show amplitude spectra of the B filter data 
with consecutive prewhitening by the low frequencies. For 
the 2000 data we calculated the window function as the 
Fourier transform of a single sinusoid with frequency 1.129 



Site 


Run start 
HJD - 2450000 


Length 
(h) 


t 

(s) 


N 


Observer 


SAAO 


1501.250 


3.02 


600 


16 


LAB 


SAAO 


1502.268 


2.39 


600 


13 


LAB 


SAAO 


1503.262 


2.39 


600 


13 


LAB 


SAAO 


1504.249 


2.99 


600 


16 


LAB 


MS 


1505.029 


0.33 


1200 


1 


DJJ 


SAAO 


1505.255 


2.42 


600 


13 


LAB 


SAAO 


1506.250 


0.35 


600 


2 


LAB 


SAAO 


1507.297 


1.29 


600 


7 


LAB 


SAAO 


1508.249 


2.81 


600 


15 


LAB 


SAAO 


1509.251 


2.54 


600 


12 


LAB 


.MS 


1509.960 


0.68 


1200 


3 


DJJ 


MS 


1510.924 


2.22 


1050 


6 


DJJ 


SAAO 


1511.256 


2.42 


600 


13 


LAB 


MS 


1511.942 


1.62 


1050 


5 


DJJ 


SAAO 


1512.254 


2.44 


600 


13 


LAB 


ESO 


1762.741 


2.70 


975 


8 


TM 


ESO 


1763.700 


3.69 


1100 


10 


TM 


SAAO 


1765.347 


0.17 


600 


1 


CK 


ESO 


1765.744 


3.69 


1100 


8 


TM 


ESO 


1766.703 


2.92 


1100 


7 


TM 


SAAO 


1768.347 


8.08 


1095 


21 


CK 


ESO 


1768.802 


1.91 


1100 


5 


TM 


SAAO 


1769.336 


8.21 


1016 


24 


CK 


SAAO 


1770.325 


6.95 


1000 


21 


CK 


SAAO 


1771.345 


8.12 


975 


24 


CK 


SAAO 


1772.345 


7.70 


1200 


20 


CK 


SAAO 


1773.354 


7.70 


1000 


23 


CK 


SAAO 


1774.406 


1.34 


1400 


3 


CK 


SAAO 


1776.328 


0.90 


1500 


2 


CK 


SAAO 


1777.388 


1.23 


1600 


3 


CK 


SAAO 


1778.348 


5.55 


1557 


7 


CK 


SAAO 


1849.267 


2.42 


600 


12 


LAB 


SAAO 


1850.271 


2.56 


600 


11 


LAB 


SAAO 


1851.274 


2.50 


600 


10 


LAB 


SAAO 


1852.262 


2.81 


600 


15 


LAB 


SAAO 


1853.258 


2.68 


600 


11 


LAB 


SAAO 


1855.354 


0.55 


700 


2 


LAB 


SAAO 


1856.258 


2.44 


600 


13 


LAB 


SAAO 


1857.265 


2.47 


600 


13 


LAB 


SAAO 


1858.264 


2.25 


600 


12 


LAB 


SAAO 


1859.263 


2.21 


600 


12 


LAB 


SAAO 


1860.265 


2.23 


600 


12 


LAB 


SAAO 


1862.262 


2.30 


600 


12 


LAB 


Total 




128.19 




476 





cycles/day and amplitude of 50 mmag. Due to the multi- 
site coverage, aliasing is not a problem. We detected six 
significant frequencies in the light variations. The first four 
frequencies were already shown in Table 1. 

The significance of a detection was estimated using the 
S/N criterion of Breger et al. (1993). The residuals after 
prewhitening this frequency solution suggest that more peri- 
odicities are present. The highest peaks in the residual spec- 
tra are located at the same frequencies in each of the B, V, 
and I c data sets, and the "noise level" decreases from B to I c . 
This suggests that most of the "noise" in the blue may well 
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Figure 4. Top panel: Spectral window of the B-filter time-series 
photometry of HD 209295 obtained in the year 2000. The other 
panels show the amplitude spectra of the data with consecutive 
prcwhitcning. Six significant periodicities are present. 



be due to additional periodicities. The derived frequencies 
and amplitudes are listed in Table 6. 



3.1.2 The 8 Scuti pulsations 

The zero point uncertainties in the data set from 1999 have 
little effect in the frequency range in which the 8 Scuti 
pulsations are present. Consequently, we can incorporate 
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Figure 5. The spectral window of the combined residual B fil- 
ter data and amplitude spectra in the frequency domain of the 
<5 Scuti pulsations. Only one frequency is convincingly detected, 
but the presence of several low-amplitude modes in a range of 12 
- 25 cycles/day is suspected; some possible further frequencies 
are indicated. 

these observations into the frequency analysis. We there- 
fore prewhitened the data by the six-frequency solution dis- 
cussed above from the 2000 data. For the 1999 data we 
prewhitened by the first five of these frequencies which were 
found in common with the analysis of the 2000 data. This 
prewhitening is important because low-frequency variations 
can artificially increase the noise level in the high-frequency 
domain through spectral leakage. The amplitude spectrum 
of the combined B filter data and subsequent prewhitening 
are shown in Fig. 5. 

Although several peaks in the amplitude spectra in 
Fig. 5 are prominent, only one signal is a significant de- 
tection. This is surprising, as typical multi-periodic beating 
is seen in the light curves of Fig. 3, e.g. HJD 2451770.35 vs. 
HJD 2451770.55. We therefore checked whether these sig- 
nals could be non-coherent by analysing several subsets of 
data. No evidence for short-lived 8 Scuti type variations or 
variations with changing amplitude was found. In fact, sev- 
eral of the more conspicuous signals were present in all data 
subsets (e.g. in different subsets of the multi-site data or 
in individual filters). Thus we suspect that several such pe- 
riodicities may be real. A multi-frequency solution with /a 
and the three next strongest variations in the 8 Scuti regime 
fits the short-term variations reasonably. The corresponding 
peaks are marked with arrows in the lowest panel of Fig. 5. 
The frequency range of possibly excited S Scuti pulsations 
can be constrained to 10 - 32 cycles/day. We list the results 
of our frequency analysis of the photometric data in Table 6. 

We need to make two remarks here. Firstly, the results 
for the data from 1999 are only listed for completeness in 
Table 6. We suspect that the amplitudes of the 7 Dor pul- 
sations in this data set are artificially decreased not only 
because of zero-point problems, but also due to poor data 
sampling. Secondly, we stress that these are only preliminary 
results. We will revisit the frequency analysis in Section 3.4. 
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Table 6. Results of a multi-frequency solution for HD 209295 derived from the photometric data. Frequency error estimates range from 
± 0.0001 cycles/day to ± 0.002 cycles/day, for the strongest and weakest 7 Dor mode respectively, and ± 0.005 cycles/day for the 5 Scuti 
mode. Signal-to-noise ratios following Brcgcr et al. (1993), calculated from the B filter data from the year 2000 are also given. A S/N > 4 
is taken to be a significant detection of a signal. 









Measurements from the year 2000 


Measurements from the year 1999 


ID 


Frequency 


S/N 


B Amplitude 


V Amplitude / Amplitude 


B Amplitude 


V Amplitude 




(d- 1 ) 




(mmag) 


(mmag) (mmag) 


(mmag) 


(mmag) 








± 0.4 


± 0.3 ± 0.2 


± 0.6 


± 0.5 


7 Doradus frequencies 


h 


1.1296 


35.9 


49.6 


38.5 23.0 


35.9 


26.9 


h 


2.3024 


28.3 


39.1 


28.3 15.5 


24.5 


17.3 


h 


2.5758 


11.6 


16.0 


11.4 6.2 


21.7 


15.8 


U 


1.1739 


9.0 


12.4 


7.9 4.6 


5.4 


4.9 


h 


1.7671 


7.0 


9.7 


7.4 4.1 


11.0 


8.5 


h 


2.2572 


5.4 


7.4 


5.3 3.1 


n/a 


n/a 


8 Scuti frequency 


/A 


25.9577 


5.3 


1.8 


1.4 0.8 


1.5 


1.3 



3.2 Spectroscopic analyses 

The first step in the analysis is to rectify the spectra, i.e. to 
place the continuum. This was done by using an unbroad- 
ened synthetic spectrum with T c ff = 7500 K, logp = 4.00 
as a template, using the SPECTRUM code (Gray & Corbally 
1994). A running median of each echelle order was divided 
by the corresponding section of the synthetic spectrum and 
taken to represent the response function of the instrument. 
A polynomial of degree 5 was fitted to the response function 
and used to correct the observed spectrum. The result is the 
rectified spectrum of the star used for cross-correlation. 



3.2.1 Radial velocities 

For each order, the observed rectified spectrum was corre- 
lated with the corresponding section of the synthetic spec- 
trum after removing the unit continuum. The resulting cor- 
relation function is, in effect, the mean line profile with the 
continuum removed. The "radial velocity" for each order is 
obtained by fitting a quadratic to the correlation function 
and finding the position of the maximum. The mean radial 
velocity from all the orders is also obtained. A standard error 
of typically 1-3 km s _1 was found for the SAAO spectra. 

A few of the spectroscopic observations listed in Table 
5 turned out to be of a quality too poor for further use. 
Four spectra were discarded, and 472 were retained. Some 
slight deviations (< 0.5 km/s) in the seasonal radial velocity 
zero-points in the SAAO spectra were rectified by means of 
telluric lines. 

An immediate finding from our time-resolved spec- 
troscopy is that the radial velocities of HD 209295 are 
strongly variable. The total radial velocity amplitude is in 
excess of 100 km/s, which cannot be due to the pulsations of 
the star. HD 209295 must be a member of a binary system. 

We therefore attempted to determine the orbital period. 
Visual inspection of the radial velocity curve suggests P or b ~ 
3 days, and Fourier analysis (upper panel of Fig. 6) indeed 
implies an orbital frequency near 0.32 cycles/day, but our 
data set is affected by aliasing. Prewhitening trial frequen- 



cies suggested that the shape of the orbital radial velocity 
curve is not sinusoidal. 

This is a situation in which use of the residualgram 
method (Martinez & Koen 1994) is indicated, which is based 
on a least-squares fit of a sine wave with M harmonics. The 
residual sum of squares at each trial frequency is evaluated; 
M can be chosen freely. In that way, alias ambiguities can be 
circumvented by taking advantage of the information in the 
Fourier harmonics. Preliminary trials with Fourier analysis 
suggested that M = 3 is a good choice for a residualgram as 
shown in the lower panel of Fig. 6. 

The residualgram method eliminated the aliasing prob- 
lem for the determination of the orbital frequency. We use 
its result for a refinement of the orbital frequency by using 
Period 98 (Sperl 1998). This leads to an improved orbital 
frequency of / orb = 0.32198 cycles/day (P orb = 3.10575 d). 

To determine the orbital solution from our radial ve- 
locities, we firstly used the SAAO data only, which are the 
most extensive (almost 90% of all spectra) and are homoge- 
neous. We weighted the individual measurements based on 
their standard errors. Measurements with standard errors 
smaller than 2 km/s were given a weight of 1; the weights 
decreased down to 0.2 for a few measurements with standard 
errors between 8-9 km/s. We determined an initial orbital 
solution from these data with an updated version (Strass- 
meier, private communication) of the differential-correction 
method (Barker, Evans & Laing 1967). 

Using the initial parameters from these methods, we ex- 
amined the relative zero-points of the three blocks of SAAO 
measurements; no third body in the system was detected. 
We then fitted all the campaign data with these orbital pa- 
rameters and examined the zero-points of the measurements 
from the other sites. The ESO radial velocities needed to be 
shifted by +7 km/s. We then calculated our final solution 
from all the data, with the mean radial velocity fixed. This 
orbital solution is shown in Table 7 and Fig. 7. 

We note that an extrapolation of our orbital solution 
to the measurements of Grenier et al. (1999) quoted in the 
introduction results in a good fit. The two measurements 
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Frequency (cycles/day) 

Figure 6. Upper panel: Fourier amplitude spectrum of the radial 
velocities of HD 209295. An orbital period near 3 days is implied, 
but severe aliasing is present. Lower panel: residual sum of squares 
spectrum of a 3-harmonic fit (RSS3) to the radial velocities. Al- 
though some power leaks into sub-harmonics, the correct orbital 
frequency is now unambiguously detected. 



Table 7. The orbital solution from radial velocities of HD 209295. 
The symbols have their usual meanings; w is the argument of 
periastron. 



Parameter 


Unit 


Value 




(days) 


3.10575 ± 0.00010 


7 


(km/s) 


-23.7 ± 0.4 


K 


(km/s) 


54.2 ± 0.7 


e 




0.352 ± 0.011 


U! 


(degrees) 


31.1 ± 2.0 


T 


HJD 


2451771.864 ± 0.014 


ai sin i 


R Q 


3.11 ± 0.04 


f(m) 


M Q 


0.042 ± 0.002 



by these authors were taken at similar orbital phases, which 
explains the comparably small radial velocity difference. 

Having solved for the orbital parameters, we can now 
discuss two issues: the nature of the companion and the pos- 
sibility of mass transfer in the system. In the upper panel of 
Fig. 8., we show the dependence of the secondary mass on 
orbital inclination for various values of the primary mass. 
From this, we see that M2 > 0.62Mq. The lower panel of 
Fig. 8 shows that HD 209295 may come close to filling its 
Roche Lobe (calculated with the approximation by Eggle- 
ton 1983) at periastron if the companion has a small mass. 



a 



o 
o 

1) 
> 



K 




0.5 1 
Orbital Phase 



Figure 7. Radial velocities of HD 209295 phased with the orbital 
solution (solid line). Open circles are SAAO radial velocities, filled 
circles are ESO radial velocities and star symbols denote the MS 
measurements. We note that although the number of spectra from 
ESO and MS is small, they cover important phases of the orbit 
and are therefore quite valuable. The formal radial velocity errors 
quoted in Sect. 3.2.1 underestimate the true errors. 



However, we will show in Sect. 4.1 that the companion is 
probably too massive for the primary to fill its Roche Lobe. 

Assuming that the companion of HD 209295 is a main- 
sequence star, its absolute magnitude M v > 4.4; otherwise 
its spectral lines would be detected. A spectral type of GO 
and later is therefore suggested with a mass < 1.05Mq 
(Drilling & Landolt 2000). Our photometry also shows no 
eclipses (within a conservative limit of 4 mmag) , allowing us 
to infer the constraint i < 80° . A main-sequence companion 
of O.62M < M 2 < 1.05Mq should be detected by means 
of its infrared excess. We will examine this in Sect. 3.5. 

After removal of the orbital variation, the residual ra- 
dial velocities can be searched for pulsational signals. The 
amplitude spectrum of these residuals can be found in Fig. 
9. Although there are signals in the same frequency range 
as in the photometry, none gives a reliable detection accord- 
ing to the signal-to-noise criterion of Breger et al. (1993). 
Consequently, we need to continue the examination of the 
spectroscopic evidence of the pulsations of HD 209295 by 
looking at the line-profile variations. 



3.2.2 Projected rotational velocity 

The correlation profiles for each season can be co-added (af- 
ter correcting with our orbital solution) to form what is es- 
sentially a mean line profile. The projected rotational veloc- 
ity, vsmi, can be determined by fitting a model profile of 
a rotating star. The model is spherical with a linear limb 
darkening coefficient u — 0.3 and the visible hemisphere is 
divided into a large number of surface elements. The in- 
trinsic line profile for each element summed over the hemi- 
sphere, weighted according to limb darkening, is computed. 
The projected rotational velocity and mean radial velocity 
are adjusted until a best fit is obtained to the observed pro- 
file. We obtain vsmi = 108 ± 10 km s _1 for the SAAO 
measurements from 1999 and the second set from 2000, and 
vsmi = 98 ± 10 km s _1 for the first data set of the year 
2000. Because of the uncertainty in the wings of the corre- 
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Figure 8. Upper panel: the mass of the companion of HD 209295 
versus orbital inclination. The solid line is for a primary mass of 
1.84 Mq as inferred in Sect. 1.1, the dotted lines represent the ltr 
limits of this determination. Lower panel: the Roche Lobe radius 
of the HD 209295 primary at periastron depending on secondary 
mass, again for M = 1.84 ± 0.07MQ. The filled circle with the 
error bar shows our radius determination for HD 209295 from 
Sect. 1.1. 
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Figure 9. The amplitude spectrum of the residual radial veloci- 
ties of HD 209295 after removing the orbital solution. 



lation profile, the fit was performed using the part of the 
profile within 90 km s _1 of the line centre. 



3.2.3 The high- order profile variations 

In a recent paper on the S Set star 38 o 1 Eri, Balona (2000) 
was able to show that pulsation modes of high spherical 
degree £ are easy to resolve in stars with moderate to high 
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Figure 10. Left panel - correlation function; right panel - grey- 
scale difference after removing the correlation function. The data 
are from the 1999 season. The abscissa is in km s and the 
ordinate is the HJD measured from HJD 2451500.000. 



rotational velocities. Bearing in mind the result from the 
previous paragraph, HD 209295 falls within this category. 
We will only use the SAAO spectra for the following analysis, 
as these are the most homogeneous and best sampled. 

Line profile variations due to modes of high degree are 
seen as moving sub-features in the correlation profiles of 
the spectra used to determine the radial velocities. As a 
particular sub-feature crosses the centre of the profile, it 
determines the radial velocity. A short while later another 
traveling sub-feature will appear and the radial velocity will 
suddenly change. These discontinuous jumps are partly re- 
sponsible for the rather high errors of the radial velocities. 

From the spacing of the traveling sub-features, we es- 
timate that the spherical degree is £ ^ 5. To enhance the 
profile variations, we removed the orbital motion and con- 
structed the mean correlation profile from the SAAO data. 
Each correlation profile was divided by the mean profile to 
construct "difference" profiles at the given times. The mov- 
ing features are most clearly seen when these difference pro- 
files are arranged in a time sequence and converted to a 
grey-scale image. An example is shown in Fig. 10. 

While it is clear from Fig. 10 that periodic high-order 
line profile variations are present, it is not possible to de- 
termine the frequencies from the radial velocities alone, as 
shown before. Clearly, the radial velocity is a poor indicator 
for modes of high degree. 

A common method to determine periodicities in line 
profile variations is to calculate the periodogram along wave- 
length bins throughout the line profile, or across the corre- 
lation profile. Most of the signal in this analysis is located at 
frequencies very close to the orbital period and its aliases, 
suggesting that the removal of the orbital variation did not 
work satisfactorily. However, we found the signatures of two 
previously detected modes in the line profiles, although they 
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Figure 11. Phase diagram for two modes of HD 209295 (lower 
panel) through the correlation profile (upper panel) , which is arbi- 
trarily centred at 6000A. The full circles denote the phases of the 
mode frequency, and the open circles that of its first harmonic. 
The zero-points in the phases are arbitrary. The dotted vertical 
lines indicate usini as determined in Sect. 3.2.2 translated into 
units of wavelength. This is also the position in the correlation 
profile where we determine the total phase change. 



Table 8. The total phase variation of the mode frequencies A</>, 
their first harmonics A0i, and the implied mode identifications 
following Telting & Schrijvers (1997). 



ID 


A<fo 


A4>! i 


H 




(n rad) 


(it rad) 




fi 


0.85 


-0.05 1 


1 


h 


1.06 


-0.77 1 


1 



cannot be detected without prior knowledge. These are the 
photometric 7 Doradus modes /1 and f%. 

This finding gives us the possibility to attempt a mode 
identification with the method proposed by Telting & Schri- 
jvers (1997). This technique utilises the phase change of 
a mode as it progresses through the line profile, and pro- 
vides some simple mode identification diagnostics depend- 
ing on the total phase change of the mode frequency and its 
first harmonic (e.g. A0o ~ n£). We calculated these phase 
changes for the two modes which we think are present and 
we show their behaviour through the correlation profile in 
Fig. 11. The mode identification implied by this analysis is 
summarized in Table 8. 

We note that Telting & Schrijvers (1997) have given 
error estimates for the mode identifications with this method 
with ±1 for I and ±2 for m. The identifications are plausible: 
fi and /2 have high photometric amplitude. 



3.3 Combining the spectroscopic and photometric 
information 

Comparing the frequencies of the 7 Doradus pulsations in 
Table 6 to the orbital frequency we see that the photometric 
mode /3 corresponds exactly to 8/ or b, and is, within the 
errors, consistent with 7/ or b- As noted in Sect. 1.1, the rela- 
tion /4 = /2 — fi is also present within the mode frequencies. 

Consequently, we revisited the frequency analysis of the 
photometric data. We again restricted ourselves to the data 
from 2000. This time, however, we used Period 98's capa- 
bility to fix signal frequencies to certain values, in our case 
integer multiples of the orbital frequency, and to perform 
nonlinear least squares fits with frequencies thus fixed. 

We started by using the six 7 Doradus frequencies in 
Table 6, and proceeded by the usual prewhitening proce- 
dure. The corresponding amplitude spectra can be found in 
Fig. 12; they are generated by combining the B filter residu- 
als with the V filter residuals multiplied by a factor of 1.29. 

The highest-amplitude peak in the upper panel of 
Fig. 12 does, indeed, correspond to 3/ or b- In addition 5/ or b 
(second panel of Fig. 12) and 9/ or b (third panel of Fig. 12) 
are also found. This latter frequency corresponds to exactly 
/1 + /s. We therefore fixed this frequency combination to its 
exact value before proceeding. After prewhitening these fre- 
quencies, some peaks, including one at 4/ orb , are still visible 
in each filter, but not at a level considered significant. 

The third panel from the bottom of Fig. 12, shows the 
residual amplitude spectra in a wider frequency range af- 
ter prewhitening by the 5 Scuti mode, /a- Again, there is 
evidence of additional 7 Doradus and weak 5 Scuti periodici- 
ties. The referee suggested to filter out the slow variability to 
examine the claimed S Scuti variability further. This is done 
in the second panel from the bottom of Fig. 12; albeit some- 
what affected by the filtering, the additional S Scuti vari- 
ations remain with similar amplitudes. They are therefore 
not artifacts from spectral leakage of the low frequencies. In 
the lowest panel of Fig. 12 we show the amplitude spectrum 
of the comparison star data. Although the number of obser- 
vations of each comparison star is half that of the program 
star, the noise level in this periodogram is lower still, again 
suggestive of further periodicities below the detection level. 

In Table 9 we show the final multi-frequency solution 
for the 2000 data set. The corresponding results for 1999 are 
of little value, as the data sampling is too poor to attempt 
such a fit. Error estimates are formal values following Mont- 
gomery & O'Donoghue (1999). While they should be quite 
realistic for the 5 Set modes, they probably underestimate 
the errors for the 7 Dor modes by a factor of 3-4. 

The rms residual errors per single data point for this 
solution are 7.0 mmag in the B filter, 5.7 mmag in the V 
filter and 4.5 mmag in the I filter. These values again suggest 
undetected periodicities, since the rms scatter of a single 
comparison star measurement is 5.1 mmag in B, 4.0 mmag 
in V and 4.2 mmag in I. 

Another estimate of the orbital period can be made 
from the frequency analysis. We find P orb = 3.1057T0.0010, 
consistent with the determination in Sect. 3.2.1. We have ex- 
amined the residual light curve for possible ellipsoidal varia- 
tions, but found none exceeding 2 mmag. Finally, one might 
suspect that f@ — 2/i , but a fit assuming this relation gives 
a significantly poorer solution than the one listed in Table 9. 
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Table 9. The final multi-frequency solution for HD 209295 from all the year 2000 photometry. Pulsational phases for mean light level 
are given with respect to a time of periastron passage, HJD 2451771.864. Amplitude signal-to-noise ratios are calculated from Fig. 12. 



ID 


Combination 


Frequency 


S/N 


B Ampl. 


B phase 


V Ampl. 


V phase 


/ Ampl. 


/ phase 






(d-l) 




(mmag) 


(°) 


(mmag) 


(°) 


(mmag) 


(°) 










± 0.3 




± 0.2 




± 0.2 




7 Doradus frequencies 


h 




1.12934 ± 0.00005 


63.7 


50.3 


-39.4 ± 0.3 


38.9 


-40.3 ± 0.3 


23.5 


-42.6 ± 0.4 


h 




2.30217 ± 0.00006 


49.9 


39.4 


-115.4 ± 0.4 


28.7 


-116.4 ± 0.4 


15.4 


-117.0 ± 0.6 


h 


8/orb 


2.57593 ± 0.00011 


23.2 


18.3 


66.5 ± 0.8 


13.2 


67.2 ± 0.9 


7.3 


64.3 ± 1.3 


h 


h-h 


1.17283 ± 0.00004 


14.4 


11.4 


-80.4 ± 1.3 


7.6 


-73.3 ± 1.5 


4.0 


-68.5 ± 2.4 


h 


9/orb — /l 


1.76859 ± 0.00005 


13.6 


10.8 


96.5 ± 1.3 


8.2 


95.4 ± 1.4 


4.8 


86.2 ± 2.0 


fo 


7/orb 


2.25394 ± 0.00011 


10.6 


8.4 


2.0 ± 1.7 


6.6 


-1.0 ± 1.8 


3.5 


2.0 ± 2.8 


fr 


3/orb 


0.96597 ± 0.00011 


8.9 


7.0 


-39.2 ± 2.0 


6.2 


-35.9 ± 1.9 


4.9 


-33.3 ± 1.9 


h 


5/orb 


1.60996 ± 0.00011 


5.8 


4.6 


-162.1 ± 3.1 


3.9 


-158.6 ± 3.0 


1.9 


-173.9 ± 5.2 


h 


9/orb 


2.89792 ± 0.00011 


5.7 


4.5 


47.0 ± 3.2 


3.5 


48.9 ± 3.3 


2.2 


53.5 ± 4.3 










8 Scuti 


frequency 










fA 




25.9577 ± 0.0015 


5.2 


1.8 


158.9 ± 8.3 


1.4 


162.8 ± 8.3 


0.8 


179.1 ± 12.2 



/orb 0.32199 ± 0.00011 



3.4 Attempts at mode identification from colour 
photometry 

The main reason for observing HD 209295 in more than one 
filter was to attempt a mode identification. This method 
relies on comparing theoretical amplitude ratios and phase 
differences of the variations in different wavebands with the 
observed values. The high amplitude of the light variations 
of HD 209295 make the star very well suited for this method. 
Because of the very low amplitude of the <5 Scuti pulsations, 
only the 7 Doradus modes are considered. However, our cal- 
culations indicate that the amplitude ratios and phase dif- 
ferences for the S Scuti mode are consistent with pulsation. 

We first attempted to apply the method developed by 
Koen et al. (1999) to the 7 Doradus modes. This technique 
uses the observed amplitude ratios and phase differences 
in all filters simultaneously. Unfortunately, no meaningful 
results were obtained. We then constructed two-colour di- 
agnostic diagrams showing amplitude ratio as a function 
of phase difference (Watson 1988). In such diagrams, the 
theoretically determined areas of interest are specified and 
compared with the observations. We used the Warsaw-New 
Jersey stellar evolution and W.A. Dziembowski's NADROT 
pulsation code for models of 1.8 and 1.9 Mq with effective 
temperatures between 7250 and 7800 K. These models span 
the possible range of parameters for HD 209295. We then 
computed theoretical amplitude ratios and phase shifts for 
the eigenmodes of those models following Balona & Evers 
(1999) and examined their location in the corresponding di- 
agrams. We noticed that (except for a few combinations of Q 
and £) the results clustered in well-defined regions depending 
on i. Hence, we defined those regions as our regions of in- 
terest and compared them to the observations. An example 
is shown in Fig. 13. 

Before evaluating the potential for mode identifications 
from Fig. 13, two comments need to be made. Firstly, we 
only show the results for the B and I c filter pair, as the 
diagrams are similar in appearance for the other two possible 
filter combinations, but the observational errors are larger. 



Secondly, the locations of the areas of interest depend 
quite strongly on the pulsation constant Q. This is a well- 
known result (e.g. discussed by Garrido 2000), and is the 
reason why we have separated the results into three subsets 
which depend on Q. We chose a suitable subdivision based 
on the observed 7 Dor frequencies. 

The implications from Fig. 13 can be summarized as 
follows. For modes with Q > 0.2 d, the mode identification 
method gives meaningless results. For the modes with Q < 
0.2 d there is some agreement, but this might be coincidence. 
In any case, the amplitude ratios of the photometric modes 
are consistent with I = 1 or 2, but not with I = 3. 



3.5 The infrared colours and the companion of 
HD 209295 

In Section 3.2.1 we have shown that the binary companion 
of HD 209295 must have a mass of at least 0.6 Mq. If it 
were a main-sequence star, this would correspond to a spec- 
tral type of K5. The most luminous possible main-sequence 
companion to HD 209295 is a GO star. A star more luminous 
than this would be detected in the spectrum. We decided to 
search for a companion in the infrared. For this purpose we 
show the system's optical and infrared colours in Table 4. 
We note that interstellar reddening is insignificant because 
of the proximity of the star (d = 122 pc) and its high galactic 
latitude (b = —44°); there is also no evidence for reddening 
in the observed Stromgren colours (Table 2). 

We adopted standard relations for absolute magnitude 
as well as optical and infrared colours of main-sequence stars 
from the tables of Tokunaga (2000) and Drilling & Landolt 
(2000). Given the observed V magnitude and (V — ic) colour 
of the system, we added the fluxes of hypothetical GO, K0 
and K5 main-sequence companions. For the primary's flux 
we used a suitable interpolation between A7 and F0 to re- 
produce the observed V and (V — i c ). The infrared colours 
of those models were then compared to observations. The 
result, displayed in Fig. 14, shows the relative infrared ex- 
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Figure 12. Top panel: the amplitude spectrum of the combined 
B and scaled V-filter residuals of HD 209295 after subtraction of 
the frequency solution of Table 6. Another harmonic of the orbital 
frequency is detected. In the following three panels we show the 
effects of further prewhitening. All new detections are again har- 
monics of / or b- Third panel from bottom: the residual amplitude 
spectrum in a wider frequency range. Second panel from bottom: 
the residual amplitude spectrum after filtering the residual low- 
frequency variations. Bottom panel: amplitude spectrum of the 
differential magnitudes of the comparison stars. 
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Figure 13. Comparison of theoretical and observed B/I c colour 
amplitude ratios and phase differences with observations. Up- 
per panel: variations with frequencies less than 1.2 cycles/day. 
Middle panel: variations with frequencies between 1.6 to 1.8 cy- 
cles/day. Lower panel: variations with frequencies greater than 
2.2 cycles/day. The areas outlined by the full lines are theoreti- 
cal predictions for £ = 1, and areas represented by dashed lines 
correspond to I = 2. Symbols with error bars are the observed val- 
ues. Filled circles correspond to independent modes, open circles 
to harmonics of the orbital period and diamonds to combination 
frequencies. Note the different abscissa scale in the middle panel. 
Regions of interest were also calculated for I = 3, but they are 
off scale towards higher amplitude ratios. The 1 = 2 area for 
0.14 < Q < 0.19 is off scale towards larger positive phase shifts. 
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Figure 14. A search for infrared excess in the energy distribution 
of HD 209295. The solid horizontal line represents an interpola- 
tion of standard relations for A- type stars to the observed (V — I c ) 
colour of HD 209295. The dotted line is the same with a GO 
main-sequence companion added, and the dashed line shows the 
expected colours if a K5 main-sequence companion were present. 
A K0 companion would be practically indistinguishable from a 
GO star in this diagram. The filled circles with the error bars are 
the actual observations. 



cess of the observations and of some of the binary models 
compared to the expected infrared magnitudes. 

Whereas the infrared J and H magnitudes seem to indi- 
cate the presence of a companion, the K magnitude argues 
against this idea; the L magnitude is too uncertain to be 
useful. The infrared excess we calculated also depends on 
the uncertainties of the standard relations and the accuracy 
of our measurements. We conclude that we cannot pinpoint 
the nature of the companion of HD 209295 from the infrared 
data alone. 



4 FURTHER DISCUSSION OF THE 
OBSERVATIONS 

4.1 The secondary component and evolutionary 
history of the HD 209295 system 

We were unable to detect the secondary component of the 
HD 209295 system in the infrared. However, we can still 
obtain further constraints on its nature. The first constraint 
is the absence of ellipsoidal variability. 

Morris (1985) derived theoretical expressions for sys- 
tem parameters of ellipsoidal variables from the observed 
light curves. Using his Eq. 6, a limb darkening coefficient 
of 0.59 from Claret (2000), a gravity darkening exponent 
of 0.84 (Claret 1999), and a conservative upper limit of 5 
mmag on the peak-to-peak amplitude of possible ellipsoidal 
variations of HD 209295, we can determine the minimum 
orbital inclination of the system as a function of secondary 
mass. The result is shown in Fig. 15. 

Combining this information with the dependence of sec- 
ondary mass on inclination from the mass function of the 
binary (upper panel of Fig. 8), we can derive a refined lower 
limit of the secondary mass. The result is M2 > 1.04Mq. 
This lower limit on the secondary mass has two implications. 
Firstly, it means that the primary is quite well within the 
Roche limit. Secondly, the secondary mass is very close to 
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Figure 15. The maximum inclination of the binary orbit of HD 
209295, consistent with the absence of ellipsoidal variability, as a 
function of secondary mass. The full line corresponds to a primary 
mass of 1.84 Mq, whereas the dotted lines delineate our lu error 
estimate of ± 0.07 Mq. 



the upper limit 1.05Mq derived in Sect. 3.2.1 assuming a 
main-sequence secondary. In addition, the orbit of the sys- 
tem is also quite eccentric, which is surprising as one would 
have expected circularization to have taken place. 

Claret, Gimenez & Cunha (1995) and Claret & Cunha 
(1997) compared theoretical orbital circularization and syn- 
chronization times for main-sequence binaries with obser- 
vations. From their results and using the surface gravity of 
HD 209295 as an age indicator, we conclude that the or- 
bital eccentricity of HD 209295 is inconsistent with normal 
main-sequence binary evolution. For a star similar in age 
to HD 209295 (log t w 8.86 yr as inferred from comparison 
with pulsational and evolutionary models), one would expect 
e < 0.1. There is reason to believe, therefore, that the sec- 
ondary star may be a degenerate object. A white dwarf may 
have a mass in excess of 1.04 Mq. However, in such a sce- 
nario orbital circularization would again have taken place. 
This is also expected because such a system must have un- 
dergone a previous mass transfer phase. 

We investigated the ultraviolet fluxes of the object as 
measured by the TD-1 satellite (Thompson et al. 1978) and 
compared it to the optical uby fluxes. For this purpose, we 
converted the measured uby magnitudes of HD 209295 to 
fluxes using the formulae of Gray (1998). We omitted the 
v band because it is dominated by the H$ line. A Kurucz 
model atmosphere with T e g — 7750 K (cf. Sect. 1.1) and 
log p=4.3 gave the best match. We then compared the pre- 
dicted UV fluxes from this fit with the TD-1 measurements 
(Fig. 16). There is quite a strong UV excess in the TD-1 
measurements at 1965 and 2365A compared to the model 
atmosphere, but little flux at 1565A. The shape of this ex- 
cess resembles an energy distribution of an object of T e s « 
15000 K. This is a suggestion that a white dwarf companion 
could be responsible. 

Consequently, we attempted to fit some white dwarf 
star model atmospheres (Koester, private communication) 
for log 5=8 and 10000 K < T eff < 16000 K to this excess. 
The V flux was normalized to the value predicted for a star 
at the distance of HD 209295. While we could reproduce 
the approximate shape of the UV excess, the resulting white 
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Figure 16. Comparison of optical uby and ultraviolet fluxes of 
HD 209295 to a T eS = 7750 K, log 3=4.3 Kurucz model atmo- 
sphere. A clear UV excess is notable. 



dwarf luminosity is far too small. A white dwarf star would 
only be able to generate 0.1 - 0.2 % of the observed UV 
excess. Hotter white dwarf models would be able to explain 
a larger fraction of the excess, but they would be inconsistent 
with the low TD-1 flux at 1565A. 

It is hard to imagine an astronomical object or a phys- 
ical process which could be made responsible for this UV 
excess and still be consistent with other observations. A 
less massive white dwarf is inconsistent with our mass con- 
straints. A subdwarf would be seen in the optical and is 
inconsistent with the primary being a Pop. I star. A late- 
type active companion would be detected in X-rays. We need 
to leave the reason for the UV excess unanswered; it may 
be possible that the TD-1 measurements are in error. Un- 
fortunately, the star was not observed by IUE (NASA-ESA 
1999). 

We also need to consider that the secondary may be a 
neutron star. In this case, the high orbital eccentricity would 
be a result of the supernova outburst, and the orbit would 
not be expected to circularize during the main sequence life 
span of the present primary (see Zahn 1977, whose work also 
suggests that the primary's rotation should not be synchro- 
nized with the orbital motion of the compact secondary). 

Comparing the orbital eccentricity and separation of the 
HD 209295 system with the X-ray binary population sim- 
ulations of Terman, Taam & Savage (1996) shows that the 
observed orbital parameters can be comfortably explained 
in this scenario. The undetected companion of HD 209295 
could therefore be a neutron star, and the system will evolve 
into an intermediate-mass X-ray binary after the present pri- 
mary has left the main sequence. The neutron star hypothe- 
sis is testable. The previous evolutionary phases should have 
left their marks in the chemical composition of HD 209295. 
However, the Stromgren metallicity index Smi = 0.016 and 
our spectra appear normal. 

If the unknown companion of HD 209295 were a neutron 
star, one would also expect that it would be detectable in 
X-rays. However, the star has not been detected by ROSAT 
(Voges et al. 1999, 2000) or Einstein (Moran et al. 1996). Fi- 
nally, the neutron star should have suffered an impulse from 
the supernova explosion, which might modify the space mo- 
tion of the binary system. We therefore calculated the galac- 



tic U, V, W velocities of HD 209295 from its 7 velocity from 
Table 7, its Hipparcos parallax (see Sect. 1.1) and proper mo- 
tion (fi a cos 8 = 26.84±0.59 mas, /xj = — 58.82±0.41 mas as 
measured by the Hipparcos satellite (ESA 1997)). We find 
U = -21 km/s, V = -12 km/s and W = +28 km/s (already 
corrected for the solar motion) which is not unusual for 
an early-F main-sequence star (cf. Gilmore & Zeilik 2000). 
The neutron star interpretation for the secondary in the HD 
209295 system therefore also has its weaknesses; we cannot 
at present identify this star's nature with certainty. 

4.2 Evidence for forced oscillations 

In Sect. 3.4 we reported the discovery that many of the pho- 
tometric frequencies are exact integer multiples of the orbital 
frequency. This raises the suspicion that they might not be 
free oscillation modes of the star, but are rather triggered 
by tidal interactions. Tidally induced nonradial oscillations 
have been searched for observationally, but the results have 
not been very convincing, with the possible exception of the 
slowly pulsating B star HD 177863 (De Cat 2001, Willems 
& Aerts 2001). As there are quite a number of frequencies 
which are exact integer multiples of the orbital frequency, 
there is reason to believe that HD 209295 is the best case 
for forced oscillations. Further support for this interpreta- 
tion needs to be sought. Some general information about 
tidal excitation is therefore useful. 

This effect has been studied by several authors theoret- 
ically, most often in connection with neutron star/main se- 
quence binaries. For instance, Kumar, Ao & Quatert (1995) 
showed that stellar p, f and low-order g-modes are not easily 
excitable through tidal effects, but that intermediate-order 
g-modes may be excited. These are, of course, exactly the 
modes in which 7 Doradus stars (and SPB stars) pulsate. 

Furthermore, tidally excited modes will have a shape 
adjusted to the gravitational potential of the exciting star. 
The tidal deformation can be decomposed into a linear su- 
perposition of the star's g- and p-mode oscillations (Press 
& Teukolsky 1977). Their dominating components are the 
I = 2, \m\ = 2 modes (e.g. see Kosovichev & Severnyj 1983). 
Unfortunately, we did not succeed in providing mode iden- 
tifications for one or more of the variations suspected to be 
forced oscillations. However, if we are really dealing with 
forced oscillations, they should be approximately in phase 
at periastron where the tidal force is at maximum. 

We have already listed pulsational phases of the differ- 
ent modes at a periastron passage in Table 9. The phase 
values are similar and seem to show some alignment of the 
supected tides. This is also the case for fs = 5/ or b assuming 
it is a I = 2, \m\ — 2 mode, as there is a 180° phase am- 
biguity. The variations are, however, not perfectly in phase, 
possibly due to nonadiabatic effects. 



4.3 The photometric and radial velocity 
amplitudes 

It is somewhat puzzling that there is little correspondence 
in amplitude between the photometric and radial velocity 
data. The two highest-amplitude 7 Doradus modes in the 
photometry are undetected in the radial velocities. While 
the range in frequency of 8 Scuti instability is in approximate 
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Table 10. A comparison of simultaneous V filter light and radial 
velocity amplitudes for 7 Doradus stars from the literature. 



Star 


2K 

(km/s) 


AV 
(mma 


g) 


2K/AV 
(km s _1 mag -1 ' 


I 


Rcf. 


7 Dor 


0.6 ± 0.2 


23 


± 


1 


26 ± 9 


3 


1 




2.6 ± 0.2 


27 


± 


1 


96 ± 8 


1 






1.2 ± 0.2 


13 


± 


1 


92 ± 8 


1 




9 Aur 


1.7 ± 1.2 


35 


± 


1 


48 ± 34 


3 


2, 3 




3.0 ± 0.8 


20 


± 


1 


148 ± 39 


3 


2, 3 




3.1 ± 1.1 


18 


± 


1 


172 ± 61 


? 


2, 3 


HR 8330 


5.2 ± 0.1 


15 


± 


1 


350 ± 20 


2 


4, 5 


HD 68192 


2.2 ± 0.2 


21.7 


± 


0.2 


101 ± 10 


? 


6 


HD 209295 


< 1 


38.9 


± 


0.2 


< 26 


1 


7 




< 1 


28.7 


± 


0.2 


< 35 


1 


7 



of both modes. We then obtain 2K/AV < 26 kms 1 mag 1 
for /1 and 2K/AV < 35 kms" 1 mag" 1 for f 2 . 

In Table 10 we show all available amplitude determi- 
nations from photometric and radial velocity measurements 
of 7 Doradus stars from the literature. We only considered 
data sets in which simultaneous radial velocity and pho- 
tometry was obtained to avoid being sensitive to amplitude 
variations. Only those cases where individual mode periods 
could be resolved are listed. 

It appears that there is an order of magnitude spread 
in the 2K/AV values for the different stars, but also for 
different modes of the same star. The limits we obtained for 
HD 209295 are at the lower end of the range to be found in 
Table 1. However, the lowest 2K/ AV values were obtained 
for modes of £ = 3, not for £ = 1 as suggested for the two 
independent modes of HD 209295. 

Aerts & Krisciunas (1996) offered an explanation for 
such diverse behaviour: the photometric variations arc 
mostly due to temperature variations, which hardly affect 
the radial velocities. In addition Aerts & Krisciunas (1996) 
suggest that, due to stellar rotation, toroidal corrections be- 
come important. These may be the cause for the large scatter 
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agreement in the photometry and radial velocities, there is 
no clear match in frequencies between the two data sets. 

Breger, Hutchins & Kuhi (1976) determined the ratios 
of velocity to light amplitudes in several S Scuti stars and 
found values between 2K/AV of 50 - 125 kms -1 mag -1 . 
Using this value for the S Scuti pulsations in HD 209295, we 
would expect the photometrically detected modes to have 
radial velocity amplitudes smaller than 0.2 km/s, and hence 
they would not be detected. However, we detected radial 
velocity variations with amplitudes of around 2 km/s in 
the S Scuti range, which should then generate photomet- 
ric V amplitudes of at least 15 mmag, but there is no trace 
of them. The only plausible explanation is that the short- 
period variations in the spectroscopy are not due to pulsa- 
tion. 

We can get limits to the radial velocity amplitude of 
the two independent 7 Doradus modes of HD 209295 by 
fitting the radial velocities with the photometric periods. 
This results in an upper limit of 1 km/s for the amplitude 
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Figure 17. Normalized linear growth rates, 77, for a model match- 
ing the effective temperature and surface gravity HD 209295. 
r\ = 1 means that a mode is driven throughout the whole model, 
whereas r) = — 1 means that a mode is damped throughout the 
whole model; modes with r\ > arc pulsationally unstable. Up- 
per panel: the whole computed frequency range. Lower panel: the 
7 Dor pulsation frequency domain. In certain frequency regions 
(depending on £), less damping is present. 



in 2K/AV in Table 10. High- resolution, high signal-to-noise 
spectra of the more rapidly rotating 7 Dor stars with simul- 
taneous photometry are necessary to test this hypothesis. 



5 THEORY 

We have carried out a stability analysis of the pulsations to 
check whether the inferred evolutionary state of HD 209295 
is consistent with its pulsational behaviour, and whether 
tidally induced oscillations are reasonable. 



5.1 Stability analysis 

We used the Warsaw-New Jersey stellar evolution and the 
NADROT pulsation codes to generate a series of models for 
HD 209295. We then investigated the pulsation modes for 
stability. The results for the model which best matches the 
parameters of HD 209295 (M = 1.8M Q , log T eS = 3.88, log 
L = 1.10, log gr=4.03) are shown in Fig. 17. 

The observed frequency domain of the 5 Scuti pulsa- 
tions of HD 209295 is very well reproduced. The unstable 
frequency range does not depend on £, which is also consis- 
tent with the observations. 7 Doradus modes are not driven 
in our models, as this currently requires a special treatment 
of convection (Guzik et al. 2000). In our computations we 
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Figure 18. Observed amplitude of the tidally induced radial ve- 
locity variations in a model for HD 209295. The full line corre- 
sponds to a rotational frequency / ro t = 1.852 d _1 and the dashed 
line to a rotational frequency / ro t = 1.854 d" 1 . The orbital period 
of HD 209295 is indicated by the dotted vertical line. 



used the standard mixing-length convection theory and we 
ignored the Lagrangian perturbation of the convective en- 
ergy flux when computing oscillations. Our attempts to re- 
produce the results by Guzik et al. (2000) concerning pul- 
sational instability were unsuccessful, and we plan to study 
this matter in the future. 

It is interesting to note that there are frequency regions 
in which damping is not as strong, and that they depend on 
I. We performed these computations up to I = 8, and found 
that the trend is maintained. The I = 1 and t — 2 frequency 
regions coincide with the ones actually observed. 

5.2 Tidal excitation 

We investigated the possibility of tidally excited oscillations 
in HD 209295 by determining the amplitude of tidally in- 
duced radial-velocity variations using the expressions de- 
rived by Willems & Aerts (2001). We calculated both the 
free and the forced oscillations of a 1.8 M® stellar model 
in the linear, adiabatic approximation. Since the amplitude 
may become quite large in the adiabatic approximation, we 
restricted ourselves to calculating the radial-velocity vari- 
ations at orbital periods for which the relative differences 
between the forcing frequencies of the dynamic tides and 
the eigenfrequencies of the free oscillation modes are not 
too small (for details see Willems & Aerts 2001). 

The resulting amplitudes of the tidally induced radial- 
velocity variations seen by an observer are displayed in 
Fig. 18 as a function of the orbital period. The compan- 
ion mass is assumed to be 1.5 Mq which corresponds to an 
orbital inclination of 31.2° and to a rotational frequency be- 
tween / r ot = 1.635 d -1 and / rot = 1.998 d" 1 . 

A high peak in the amplitudes is seen to occur near 



the observed orbital period of HD 209295. The peak results 
from the resonant excitation of the modes gf 4 and g^~ 6 by 
the dynamic tides associated with the forcing frequencies 
3 /orb and 4/ or b, respectively. Both tides are of spherical 
degree 1 = 2 and azimuthal order m = —2. This result 
is quite encouraging, as a signal at 3 / or b was detected in 
the photometry, and one at 4 / or b was suspected. On the 
other hand, not all of the detected orbital harmonics are 
excited in the model, but they can be reconciled with linear 
combinations of these two signals. 

We need to make several remarks here. Firstly, the 
above "solution" is probably not unique. In addition, the 
calculations have been performed using an adiabatic code, 
so they show the presence of the resonances, but not the 
real amplitudes. Secondly, attempts with lower companion 
masses (and thus lower rotational frequencies) were not very 
successful. This is interesting, as it may help to constrain the 
nature of the orbital companion of HD 209295. 

We conclude that resonant excitation is possible, but 
for precise results there are still too many degrees of free- 
dom. The prospect of developing tides as a probing tool is 
a promise for future studies. 



6 CONCLUSIONS AND OUTLOOK 

By means of 128 h of time-resolved high-resolution spectro- 
scopic and 280 h of time-series photometric observations we 
have shown that HD 209295 is not only the first pulsating 
star to be a member of two classes (a 7 Doradus as well as a 
8 Scuti star), but we have also discovered good evidence for 
the presence of tidally excited pulsation modes. The binary 
companion is probably a neutron star or a white dwarf. 

These findings make HD 209295 a rather interesting as- 
trophysical laboratory. We were able to confront photomet- 
ric and spectroscopic mode identification methods. Whereas 
the method by Telting & Schrijvers (1997) yielded reason- 
able assignments for the two independent modes dominat- 
ing the photometry (both are I = 1, |m| = 1), analysing the 
photometric colour amplitude ratios and phase differences 
did not yield any meaningful results. The necessity of using 
as many mode identification methods as possible to secure 
reliable results is stressed. 

We showed that the presence of tidally induced modes 
can be explained by theoretical models, but not all the ob- 
servations can be understood. It would be important to cal- 
culate more realistic models of tidal interaction. 

We do not fully understand the wealth of information 
provided by the variations in HD 209295. Large coordinated 
efforts will be required for further progress as our data are 
already fairly extensive. Another multi-site photometric and 
spectroscopic campaign on an even larger scale would be de- 
sirable. This should be carried over a long time base in order 
to resolve the 7 Dor pulsations, and a large data set is needed 
to detect more S Scuti modes. Large telescopes are required 
for the spectroscopy to obtain high time resolution and a 
high signal to noise ratio. This is necessary for identifying 
further modes and for detailed line-profile analysis. Photo- 
metric observations in more filters than previously obtained 
are necessary, perhaps extending the wavelength range into 
the near infrared. UV and infrared spectroscopy would be 
useful in a search for the companion star, and an abundance 
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analysis may reveal an unusual chemical composition caused 
by the possible previous evolution of the system through 
mass transfer, common envelope and/or supernova stages. 
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